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Abstract. Over the past years PHOBOS has continued to analyze the large datasets 
obtained from the first five runs of the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory. The two main analysis streams have been pursued. 
The first one aims to obtain a broad and systematic survey of global properties of 
particle production in heavy ion collisions. The second class includes the study of 
fluctuations and correlations in particle production. Both type of studies have been 
performed for a variety of the collision systems, covering a wide range in collision 
energy and centrality. The uniquely large angular coverage of the PHOBOS detector 
and its ability to measure charged particles down to very low transverse momentum 
is exploited. The latest physics results from PHOBOS, as presented at Quark Matter 
2008 Conference, are contained in this report. 

I For the full list of PHOBOS authors and acknowledgments, see appendix "Collaborations" 
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PACS numbers: 25.75.-q 

1. Introduction 

PHOBOS [1] is a heavy ion experiment at the BNL RHIC colhder, known for its 
capabihty of measuring charged particles over a broad angular acceptance, by far 
the largest of all RHIC experiments. With the PHOBOS multiplicity array charged 
particles are measured in pseudorapidity range | ^ |< 5.4 and over 27r in the azimuth. 
In addition, the unique design of the two-arm spectrometer allows for extending the 
particle momentum measurements to the lowest limit reachable at RHIC. With this 
detector we have collected data on p+p, d+Au, Cu+Cu and Au+Au collisions during 
the RHIC 2000-2005 runs, covering a wide range of collision energy and centrality. This 
comprehensive data set allowed for a systematic studies of the overall features of particle 
production mechanisms in nuclear and elementary collisions, which constituted the main 
part of the baseline PHOBOS physics program. While this part is now nearly completed, 
our current effort is mainly focused on the study of fluctuations and correlations in 
particle production, the study of which can provide deeper insight into different stages 
of the system evolution. 

In this paper we briefly summarize the results on global properties of charged 
particle production, including antiparticle to particle ratios and particle yields at 
very low transverse momentum (pt)- The emphasis is put on our recent results 
obtained from fluctuation and correlation studies. The dynamical fluctuations of the 
elliptic flow, corrected for the non-flow effects are presented for Au-fAu collisions at 
y/SNN = 200 GeV. The same high-statistics data set is used to investigate structures in 
the near- and away-side correlations with respect to high-p^ trigger particle over a broad 
range in Arj. Finally, the results from a systematic study of the two-particle angular 
correlations in p+p, Cu+Cu and Au+Au collisions at ^/sNN = 200 GeV are shown. 

2. System-size dependence of particle production 

Recently PHOBOS has completed a systematic study of the bulk properties of the 
produced particles, hke total multiplicity, dNch/drj and dN^h/dpT distributions, particle 
composition and collective flow effects, in Cu+Cu and Au+Au coUisions as a function of 
the collision energy and centrality. The results, in comparison to d+Au and elementary 
p+p collisions, show that particle production in heavy ion collisions can be described 
in terms of simple scaling rules. Number of participant (Npart) scaling is observed for 
the total particle multiplicity [2], net-proton yields [3] and also, as will be shown later, 
the low-pT particle yields. Energy and centrality dependencies of mid-rapidity yields 
factorize over an extended range of transverse momenta |4l |5]. Particle yields and 
directed (f i) and elliptic flow (^2) signals measured over a wide range of high to almost 
central pseudorapidities show energy-independence when viewed in the rest frame of 
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one of the colliding nuclei [2]. The antiparticle to particle ratios show at most weak 
dependence on the system size [6], as is illustrated by the compilation of PHOBOS data 
presented in Figure [H 
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Figure 1: Antiparticle to particle ratios for pions, kaons and protons as a function of collision 
centrality for different collision systems, measured at mid-rapidity region of 0.2 < r/ < 1.4. 

It is also observed that Npart unifies the measurements of the nuclear modification 
factors in Cu+Cu and Au+Au collisions [7]. A consistent and unified description of the 
elliptic flow measured in Cu+Cu and Au+Au collisions can be obtained after scaling w 2 
with the participant eccentricity, epart [S]- 

The above observations show that the collision geometry has a major impact on 
the dynamical evolution of the system. They also provide a tool for extrapolating RHIC 
data to the LHC energy regime, which is of particular importance since to date no theory 
or model can consistently explain the scaling rules controlling the particle production 
in heavy ion collisions. 



3. Particle production at very low transverse momenta 

The PHOBOS spectrometer design provides a unique capability of studying particle 
production at the lowest transverse momenta accessible at RHIC. The new results on 
low-pT yields for pions, kaons and protons measured in a high-statistics sample Au+Au 
collisions at ^/sNN = 200 GeV have been shown at this Conference [9]. Figure [2] shows 
the invariant low-pT yields measured for 6% of the most central Au+Au collisions in 
comparison to the PHENIX data at higher pt |iQj. The low-p-r yields agree with 
the extrapolations of Blast- Wave and Bose-Einstein fits to the PHENIX results. A 
flattening of the shape of the pt spectra is observed, stronger for heavier particles, 
which is consistent with a rapid transverse expansion of the system. Invariant yields, 
integrated over the low-p-r range (0.031 — 0.053 GeV/c for pions, 0.105 — 0.128 GeV/c 
for kaons and 0.143 — 0.206 GeV/c for protons and antiprotons) and normalized per 
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Npart, are shown as a function of Npart for Au+Au collisions at 200 and 62.4 GeV [3] 
in Figure [3l One can see that, within the errors, the invariant low-p^ yields scale with 
Npart- For more details see Ref.|9]. 
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Figure 2: Identified particle px spectra Figure 3: The integrated low-px yields for 
near mid-rapidity in Au+Au collisions at Au+Au collisions at 200 GeV (solid symbols) 

and 62.4 GeV (open symbols) as a function of 

centrality. 



^/s^ = 200 GeV. 



4. Elliptic flow fluctuations 



The measured magnitude of the elliptic flow signal, V2, should reflect the initial spatial 
eccentricity of the overlap region of the colliding nuclei. As PHOBOS has shown 
[H], the participant eccentricity, epart, which takes into account the fluctuations in the 
participant nucleon positions in the overlap region, is the relevant eccentricity which 
drives the observed azimuthal anisotropy. We have used the Monte Carlo Glauber 
(MCG) approach to study the robustness of the participant eccentricity and also to 
calculate the higher order cumulants of epart quantifying the event-by-event fluctuations 
of the initial source eccentricity [H]. It was shown that the participant eccentricity is 
a robust quantity, insensitive to the choice of MCG parameters and model assumptions 
and that it strongly fluctuates with a{epart)/ impart) varying from about 35% for peripheral 
up to 50% for most central Au+Au collisions. 

Comparison of the V2 measurements to hydrodynamic model predictions [12] 
indicate that the matter created in heavy ion collisions at RHIC has properties 
resembling those of a perfect fluid [13j. In hydrodynamic models, the fluctuations in 
the shape of the initial collision geometry should be reflected in the fluctuations of V2- 
PHOBOS has measured large (of the order of 40 — 50%) relative fluctuations of the 
elliptic flow, cr(i;2)/(^2) for Au+Au collisions at ^/sNN = 200 GeV [T3]. The main 
difficulty in this measurement is to disentangle dynamical fluctuations of f 2 and non- 
flow correlations due, for example, to resonance decays or jet production. PHOBOS has 
developed a new method to extract the non-flow component from the data. The method 
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relies on the determination of the second Fourier coefficient of two-particle angular 
correlations measured over a broad range of Ar^, thanks to the large acceptance of the 
PHOBOS detector. This coefficient, V2{r]i,T]2), contains contributions from genuine flow 
correlations, f2('7i) x V2{r]2), and non-flow correlations, S{r]i,r]2). At Arj > 2, we assume 
that the non-flow component is small. In fact a small non-flow effect at Arj > 2 was 
estimated from HIJING simulations and accounted for. Thus, performing the fit to the 
second Fourier coefficient of two-particle angular correlations allows us to determine 
the flow component, V2{f]i) x ^2(772), and then estimate the non-flow contribution by 
subtracting flow component from f|(?7i, 772). The contribution from non-flow correlations 
to cr(t>2) is J{5)/2 [15]. Figure m shows the measured total relative V2 fluctuations and 
the non-flow contribution to these fluctuations. The non-flow contribution extracted 
from the data is large, of the order of 25 — 30%. The relative flow fluctuations, after 
subtracting the non-flow contribution, are shown in Figure and compared to the 
relative fluctuations of the initial eccentricity calculated from MCG [11] and Color Glass 
Condensate (CGC) [IB] models. It can be seen that the magnitude of the relative flow 
fluctuations agree, within the errors, with both MCG and CGC calculations of the 
fluctuations in initial source eccentricity, leaving essentially no room for other, later- 
stage contributions. For more details see Ref. [T7] . 
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Figure 4: Total measured relative flow fluc- 
tuations (circles) and the contribution from 
non-flow correlations (squares) for Au+Au col- 
lisions at y^ijviv = 200 GeV. Bands show 90 % 
C.L. systematic errors. 




Figure 5: Relative flow fluctuations, corrected 
for the non-flow effects, compared to the 
relative fluctuations of the initial eccentricity 
calculated from MC Glauber and CGC 
models. 



5. High-pT triggered two-particle correlations 



Measurements of correlations with respect to the high-p^ trigger particle allow us to 
study the medium response to energetic partons produced in early hard scattering 
processes and then propagated through the dense medium created in heavy ion collisions. 
Here we report the ffist measurements of the A?7 and A0 correlations between the 
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triggered charged particle with px > 2.5 GeV/c and all associated hadrons recorded 
in the PHOBOS Octagon detector p]|. The analysis utilizes the extensive coverage in 
pseudorapidity (| rj \< 3) and azimuthal angle (A0 = 27r) of the Octagon allowing for 
the study of both short- and long-range correlations in At]. A particular emphasis is 
put on testing the presence of long-range correlations at near- and away-sides over the 
uniquely broad acceptance in Arj (—4 < Arj < 2). The correlated yields per trigger 
particle, with the elliptic flow effects subtracted, are calculated for Au+Au collisions at 
^/sNN = 200 GeV and compared to the correlated particle production in p+p events 



simulated with PYTHIA. The two dimensional correlated yield. 
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Au+Au collisions exhibits a much broader away-side peak in A(f) as compared to p+p 
correlations. The near-side yield in central Au+Au collisions shows a clear jet-like 
peak near A0 ?» as well as pronounced correlations extending over a broad At] range 
(the so called ridge structure). In Figure E] the correlated near-side yield (integrated 
over I A0 |< 1) is plotted as a function of At] for 10% of the most central Au+Au 
collisions. A prominent ridge effect is clearly visible up to At] = —4 in contrast to p+p 
PYTHIA events, where the structure is absent. The momentum kick model [18] is able 
to quantitatively describe this structure. Centrality dependence of the near-side yield, 
averaged over —4 < Arj < —2, is shown in Figure [71 One can see that the ridge effect 
weakens with decreasing collision centrality and almost completely disappears for Npart 
less than about 80. For more details of this analysis see Ref. J19j . 
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Figure 6: The correlated near-side yield for Figure 7: The near-side yield averaged over 
0-10% central Au+Au collisions compared to —4 < Ar] < —2 as a function of Npart- 
p+p PYTHIA simulations (dashed line) and 
momentum kick model predictions (solid line) . 



6. Two-particle angular correlations 



The study of the correlations between the final-state particles should provide additional 
information on the underlying mechanism of particle production at freeze-out. PHOBOS 
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has performed systematic studies of two-particle angular correlations in (A?7, A0) in 
elementary p+p collisions [20] as well as Cu+Cu collisions [21] . Recently these studies 
have been also carried out for Au+Au collisions at ^/snn = 200 GeV. The broad 
coverage of the PHOBOS Octagon detector is fully utilized in these studies. The 
multiphcity-independent two-particle correlation function, R{Ari,A(j)), is calculated 
using the method detailed in pO]. We have concentrated on studying the short- 
range correlations by projecting the two-dimensional correlation function into a one- 
dimensional correlation function R{Ari). R{Ari) exhibits clear short-range correlations 
which can be described using a simple cluster model. The one-dimensional correlation 
function is fitted to a parameterization derived in an independent cluster model [22] 
in order to extract the effective cluster size, K^/f. Figure [8] shows the effective cluster 
size as a function of Npart for ^/s^n = 200 GeV Cu+Cu and Au+Au collisions. For 
both collision systems the cluster size decreases with increasing centrality. Furthermore, 
for the same Npart, clusters in Au+Au collisions are larger than in the Cu+Cu system. 
Interestingly, this dependence on the size of the colliding nuclei disappears after plotting 
the same data as a function of fractional cross-section as shown in Figure [91 More 
differential studies of the effective size of near- and away-side clusters have also been 
performed. It is observed that away-side clusters are smaller and decrease more rapidly 
with increasing centrality than near-side clusters. This observation can be qualitatively 
explained by assuming that near-side clusters are preferentially produced close to the 
surface of the emission zone, while away-side clusters seem to propagate through the 
medium. For more details on the two-particle correlation study see [23]. 
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7. Summary and outlook 

PHOBOS is continuing to provide interesting and unique data on particle production in 
different collision systems in the RHIC energy range. The uniqueness of the PHOBOS 
data is owed to the distinct features of the PHOBOS detector: its extensive acceptance 
for charged particle measurements and capability of measuring particles at very low 
Pt- Comprehensive results on global features of the particle production are only 
briefly summarized in this report. The measurements of low-p^ yields, for the first 
time performed as a function of centrality, are presented for Au+Au collisions at 
^/sNN = 200 GeV. No anomalous enhancement of pion production is observed at very 
low pt- The spectral shapes for heavier particles show fiattening effects, consistent with 
a strong transverse expansion of the system. 

In this report we have presented a wealth of results on the fluctuations and 
correlation measurements from PHOBOS. New results on event-by-event elliptic flow 
fluctuations, corrected for the non-flow effects extracted from data, are shown for 
Au+Au collisions at y^^NN = 200 GeV. Although non-flow correlations contribute 
signiflcantly to the measured flow fluctuations, the corrected relative flow fluctuations 
are large, with a magnitude in agreement with calculations of fluctuations in the 
participant eccentricity. These results indicate that system thermalizes very rapidly 
and the initial-state event-by-event source fluctuations are efficiently converted into 
final-state momentum fiuctuations. 

The PHOBOS measurements of correlations between the trigger high-p^ particle 
and associated particles with no px cut imposed, allow for studying the ridge structure 
in the near-side correlations over the longitudinal extent larger than accessible in other 
RHIC experiments. The results for Au+Au collisions at ^/snn = 200 GeV show that the 
ridge structure of the near-side correlations persists up to the limit of our acceptance, 
i.e. Arj = 4. With decreasing collision centrality the structure becomes less pronounced, 
eventually disappearing at Npart of about 80. These results provide valuable constraints 
on models aimed to describe jet propagation through a dense medium. 

The systematic studies of two-particle angular correlations for p+p, Cu+Cu and 
Au+Au collisions at ^/sNN = 200 GeV show that at freeze-out particles tend to be 
produced in clusters with a non-trivial centrality dependence observed in nucleus- 
nucleus collisions. Furthermore, an unexpected system size dependence is seen. These 
observations provide a challenge to models describing particle production in heavy ion 
collisions. 

The PHOBOS will continue the analysis of the existing comprehensive dataset. 
The studies of the \ow-pT particle production as well as of elliptic fiow fiuctuations 
and triggered correlations will be extended to other systems and energies. Although 
PHOBOS has accomplished its initial goal, the future studies should still provide results 
improving our understanding of the physics of heavy ion collisions. 

This work was partially supported by U.S. DOE grants DE-AC02-98CH10886, DE- 
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DE-AC02-06CH11357, by U.S. NSF grants 9603486, 0072204, and 0245011, by Polish MNiSW 
grant N N202 282234 (2008-2010), by NSC of Taiwan Contract NSC 89-21 12-M-008-024, and 
by Hungarian OTKA grant (F 049823). 
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